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1 Introduction 

Factorization of systems of differential equations was first studied for the case 
of a single linear ordinary differential equation (LODE) with linear ordinary 
differential operator (LODO) of the form 

L = f (x)D n + f 1 (x)D n - 1 + ... + f n (x), D = d/dx, (1) 

where the coefficients f s (x) belong to some differential field K. Factorization 
is a useful tool for computing a closed form solution of the corresponding 
linear ordinary differential equation Ly = as well as determining its Galois 
group (see for example [IHJ HH [M])- For simplicity and without loss of 
generality we suppose that operators ([I]) are reduced (i.e. fo(x) = 1) unless 
we explicitly state the reverse. The most popular case of the differential field 
K = Q(x) of rational functions with rational or algebraic number coefficients 
is a nontrivial example which is well investigated and will be considered 
hereafter when we discuss any constructive results. 

In this paper we give a review of the current state of the theory of factor- 
ization of ordinary and partial differential operators and even more generally, 
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of systems of linear differential equations of arbitrary type (determined as 
well as overdetermined) . We start with elementary algebraic theory of factor- 
ization of linear ordinary differential operators (CQ) developed in the period 
1880-1930. After exposing these classical results we sketch more sophisti- 
cated algorithmic approaches developed in the last 20 years. The revival of 
this theory in the last two decades is motivated by the development of pow- 
erful computer algebra systems and implementation of nontrivial algebraic 
and differential algorithms such as factorization of polynomials and indefinite 
integration of elementary functions. 

The main part of this paper will be devoted to modern generalizations of 
the factorization theory to the most general case of systems of linear partial 
differential equations and their relation with explicit solvability of nonlinear 
partial differential equations based on some constructions from the ring the- 
ory, theory of partially ordered sets (lattices) and that of abelian categories. 
Many of the results of this paper may be exposed within the framework of 
the Picard-Vessiot theory. But we follow a much simpler algebraic approach 
in order to facilitate the aforementioned generalizations. 

The proper theoretical background for the simplest case — factorization of 
linear ordinary differential operators with rational coefficients — was known 
already in the end of the XIX century. Paradoxically, mathematicians of that 
epoch had developed even a nontrivial (theoretical) algorithm of factorization 
of such operators [I]! A review of this theory can be found in [41] . In 
contrast to the well-known property of uniqueness of factorization of usual 
commutative polynomials into irreducible factors, a simple example D 2 = 
D ■ D = (D + l/(x — c)) • (D — l/(x — c)) shows that some LODO may have 
essentially different factorizations with factors depending on some arbitrary 
parameters. Fortunately according to the results by E. Landau [2H] and 
A. Loewy [301 EI] exposed below all possible factorizations of a given operator 
L over a fixed differential field have the same number of factors in different 

expansions L — L\ = L\ L r into irreducible factors and the 

factors L s , L p are pairwise "similar" . (Hereafter we always suppose the order 
of factors to be greater than 0: ord(Lj) > 0, ord(Lj) > 0). We outline the 
main ideas of this classical theory in Section [2J For simplicity we discuss 
here only the case of differential operators, a generalization for the case of a 
general Ore ring (including difference and g-difference operators, see [6], E]) 
is straightforward. 

Subsequent Sections are devoted to different aspects of the theory of 
factorization of linear partial differential operators. 
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2 Factorization of LODO 



The basics of the algebraic theory of factorization of LODO was essentially 
given already in (3UJ [31] , [33] [36]. Algebraically the main results are just an 
easy consequence of the fact that the ring K[D] of LODO with coefficients 
in a given differential field K is Euclidean: for any LODO L, M there exist 
unique LODO Q, R, Qi, Ri such that 

L = Q ■ M + R, L = M ■Q 1 + R u ord(R) < ord(M), ord(Rx) < ord(M). 

For any two LODO L and M using the right or left Euclidean algorithm one 
can determine their right greatest common divisor rGCD(L, M) = G, i.e. 
L = L\ ■ G, M = Mi ■ G (the order of G is maximal) and their right least 
common multiple rLCM(L, M) = K, i.e. K = M-L = L-M (the order of K 
is minimal) as well as their left analogues IGCD and ILCM. All left and right 
ideals of this ring are principal and all two-sided ideals are trivial. Operator 
equations 

X ■ L + Y ■ M = B, L- Z + M -T = C (2) 

with unknown operators X, Y, Z, T are solvable iff rGCD(L, M) divides B 
on the right and IGCD(L, M) divides C on the left. We say that an operator 
L is (right) transformed into L\ by an operator (not necessary reduced) B, and 

write L L x , if rGCD(L, B) = 1 and K = rLCM(L, B) = L x ■ B = B x ■ L. 
In this case any solution of Ly = is mapped by B into a solution By 
of Liy = 0. Using (2) one may find with rational algebraic operations an 

operator B\ such that L\ — ^ L, B\ ■ B — l(modL). Operators L, L\ will 
be also called similar or of the same kind (in the given differential field K). 
So for similar operators the problem of solution of the corresponding LODE 
Ly = 0, Liy = are equivalent. One can define also the notion of left- 
hand transformation of L by B into L\. K = ILCM(L, B) — B ■ L\ — 
L ■ B\. Obviously left- and right-hand transformations are connected via the 
adjoint operation. Also one may prove ([35]) that two operators are left-hand 
similar iff they are right-hand similar. A (reduced) LODO is called prime 
or irreducible (in the given differential field K) if it has no nontrivial factors 
aside from itself and 1. Every LODO similar to a prime LODO is also prime. 
Two (prime for simplicity) LODO P and Q are called interchangeable in 
the product P ■ Q and this product will be called interchangeable as well if 
P ■ Q — Qi ■ Pi, Qi ^ P, Pi ^ Q. In this case P is similar to Pi, Q is similar 

to Qi and Pi P. 

Theorem 1 (Landau |28j and Loewy [30J) Any two different decompo- 
sitions of a given LODO L into products of prime LODO L — P\ Pk = 
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Pi P p have the same number of factors (k = p) and the factors are 

similar in pairs (in some transposed order). One decomposition may be ob- 
tained from the other through a finite sequence of interchanges of contiguous 
factors (in the pairs Pi ■ Pi+i). 

All definitions here are constructive over the differential field of rational 
functions K = Q(x): either using the Euclidean algorithm or finding rational 
solutions of LODO one can determine for example if two given LODO are 
similar or find all possible (parametric) factorizations of a given LODO with 
rational functional coefficients [U El [7J 150] . 

Landau-Loewy theorem also has a useful for the following ring-theoretic 
interpretation. Namely, every L e generates the corresponding left 

ideal \L); L\ divides L on the right iff \L) C \Li). If we have a factorization 
L = L\ ■ ■ ■ L k then we have a chain of ascending left principal ideals \L) C 
\L 2 ---L k ) C \L 3 ---L k ) C ... C \L k ) C 1 1) = K[D]. If the factors L k 
are irreducible, the chain is maximal, i.e. it is not possible to insert some 
intermediate ideals between its two adjacent elements. The Landau-Loewy 
theorem is nothing but the Jordan-Holder -Dedekind chain condition: 

Theorem 2 Any two finite maximal ascending chains of left principal ideals 
in the ring K[D] of LODO have equal length. 

Similarity of irreducible factors can be also interpreted in this approach. Even 
more general lattice-theoretic interpretation turned out to be fruitful for a 
generalization of this simple algebraic theory for the case of factorizations 
of partial differential operators [51]. Namely, let us consider the set of (left) 
ideals in K[Z?] as a partially ordered by inclusion set M. (called a poset). 
This poset has the following two fundamental properties: 

Property I) for any two elements A, B e M. (left ideals!) one can find a 
unique C = sup(A, B), i.e. such C that C > A, C > B, and C is "minimal 
possible". Analogously there exist a unique D = m£(A, B), D < A, D < B, 
D is "maximal possible". 

Such posets are called lattices [22j. sup(A, B) and inf(A, B) correspond 
to the GCD and the LCM in K[D}. 

For simplicity (and following the established tradition) sup(A,B) will be 
hereafter denoted as A + B and inf (A, B) as A ■ B; 

Property II) For any three A, B, C G M. the following modular identity 
holds: 

(A-C + B)-C = A-C + B-C 
Such lattices are called modular lattices or Dedekind structures. 
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As on can prove, modularity implies the Jordan-Holder-Dedekind chain 
condition: any two finite maximal chains L > L\ > ■ ■ ■ > L& > and 
L > Mi > ■ ■ ■ > M r > for a given L G M. have equal lengths: k = r 
(the same for ascending chains). For the interpretation of the notions of 
similarity, direct sums, Kurosh & Ore theorems on direct sums cf. [51 j . 

But even more fruitful for generalizations is the followings categorical 
interpretation of similarity of LODO and the Jordan-Holder-Dedekind chain 
condition. Namely, let us consider the following abelian category COVO 
of LODO. 

Objects of COVO are reduced operators L = D n +ai(x)D n ~ 1 + . . .+a n (x), 
<2j G K. One may ideally think of (finite-dimensional!) the solution spaces 
Sol(L) of such operators in some sufficiently large Picard-Vessiot extension of 
the coefficient field as another way of representation of an object in COVO. 
This helps to understand the meaning of some definitions below, but one 
shall remember that these solution spaces are not constructive unlike the 
objects-operators. 

Morphisms Hom(L, Li) in this category are constructively defined as a 
not necessary reduced LODO B such that L\ • B = C • L for some other 
LODO C. Non-constructively this B may be seen as a mapping of solutions 
of L into solutions of L\. Note that all operators here have coefficients in 
some fixed differential field K. Two operators Bi, B2 generate the same 
morphism iff B\ = B 2 (mod L). Also we should remark that this definition 

is not equivalent to the definition of a transformation of operators L — > L\ 
introduced earlier, because for morphisms: 

1) B and L may have common solutions, i.e. a nontrivial rGCD(B, L). 
This means that the mapping of the solution space Sol(L) by B may have a 
kernel Sol(rGCD(B , L)) . The morphism is not injective in this case. 

2) The image of the solution space Sol(L) may be smaller than Sol(Li). 
The morphism is not surjective in this case. 

Algebraically this means that L x ■ B = C ■ L ^ rLCM(B, P). 

Similarity of operators L and Li now simply means isomorphism of the 
objects L and L\ in this category. 

The following fact is a direct corollary of our representation of this cate- 
gory as a subcategory of the category of finite-dimensional vector spaces and 
linear mappings preserving direct sums, products etc.: 

Theorem 3 The category COVO is abelian. 

Among the many useful results for abelian categories (cf. for example [TUlfTU] ) 
we need the following 
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Theorem 4 Any abelian category with finite ascending chains satisfies the 
Jordan- Holder property. 

This will serve us as a maximal theoretical framework for an algebraic inter- 
pretation and generalization of the Landau-Loewy theorem for the case of 
systems of linear partial differential equations below. Again there exist no- 
tions of direct sums, Kurosh-Ore theorems on direct sums and the powerful 
technique of modern homological algebra for abelian categories [T6l [19]. We 
will see below that this rather high level of abstraction allows a very natural 
generalization of the definition of factorization for arbitrary systems of linear 
partial differential equations (LPDE). 



3 Factorization of LPDO 

In contrast to the case of ordinary operators, two main results are seemingly 
lost for LPDO: Landau-Loewy theorem and the possibility to use some known 
solution for factorization of operators. In the case of a LODO L obviously if 
one has its solution Lcf) = then one can split off a first-order right factor: 

L = M ■ — ^-V For a LPDO, even if one knows the complete set of 

solutions, the operator may not be factorizable, as the following classical 
examples shows: 

Example 1. The equation Lu = ^D x D y — ^ y yt j u = with D x = d/dx, 
D y = d/dy has the following complete solution: 

x + y 

where F(x) and G(y) are two arbitrary functions of one variable each. On 
the other hand, as an easy calculation shows, the operator L can not be 
represented as a product of two first-order operators (over any differential 
extension of the given coefficient field K = Q(x,y)). 

Example 2. The equation Lu = {^D x D y — ^ x + y y j u = again has the 
following complete solution: 

12(F(x) + G(y)) Q(F'(x)+G>(y)) 
U= (x + y y xTy- + F ^+G(y), (3) 

but the operator is again "naively irreducible" . More generally, the equation 

c 

Lu = u xv — — u = 0, c = const (4) 

(x + y) 2 
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has a complete solution in an "explicit" form similar to ([3]) iff c = n(n + 1) 
for n e N. In this case it has the complete solution in the form 

u = c F + Cl F' + ... + c n F^ + d G + d x G' + ... + d n+1 G^ (5) 

with some definite Ci(x,y), di(x,y) and two arbitrary functions F(x), G(y). 
Only for c = the corresponding operator is "naively reducible": L = 

D X • Dy. 

The solution technology used here is very old and can be found in [9j [20] 
under the name of Laplace transformations or Laplace cascade method: after 
a series of transformations of (jlj) one gets a naively factorizable LPDE! 

Another unpleasant example is also ascribed in [5] to E. Landau: if 

P = D x + xD y , Q = D X + 1, 
R = Dl + xD x D y + D x + (2 + x)D y , {0) 

then L = QoQoP = RoQ. On the other hand the operator R is absolutely 
irreducible, i.e. one can not factor it into product of first-order operators 
with coefficients in any extension of Q(x, y). So there seems to be no hope 
for an analogue of Landau-Loewy theorem for decomposition of LPDO into 
product of lower-order LPDO. 

We see that a "naive" definition of a factorization of a LPDO as its 
representation as a product (composition) of lower-order LPDOs lacks some 
fundamental properties established in the previous Section for factorization 
of LODO. 

Recently [51] an attempt to give a "good" definition of generalized fac- 
torization was undertaken. In the next subsection we only briefly sketch the 
ideas of this approach and describe its nontrivial relation to explicit integra- 
bility of nonlinear partial differential equations. 

3.1 General theory of factorization of an arbitrary sin- 
gle LPDO, ring-theoretic approach 

Our goal in [51] was to define a notion of factorization with "good" properties: 

• Every LPDO L shall have only finite chains of ascending generalized 
factors. In particular D x should be irreducible. 

• Jordan-Holder property: all possible generalized factorizations of a 
given operator L have the same number of "factors" in different ex- 
pansions into irreducible factors and the "factors" should be pairwise 
"similar" in such expansions. 
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• Existence of large classes of solutions should be related to factorization. 

• Classical theory of integration of LPDO using the Laplace cascade 
method should be an integral part of this generalized definition. 

An obvious extension of the definition may be suggested if one will use as- 
cending chains of arbitrary (not necessary principal) left ideals starting from 
the left ideal generated by the given operator: 

|L)c/iC/ 2 c...c/ fc c|l>. (7) 

Unfortunately one can easily see that even for the operator D x we have such 
chains, and they have unlimited length: \D X ) C \D X , D™) C \D X , D™^ 1 ) C 
. . . \D X , D y ) C |1)! So we shall take some special class of ideals, more general, 
than the principal ideals, but much less rich then arbitrary left ideals. In [51] 
we gave a definition of such a suitable subclass of left ideals called divisor 
ideals. 

For such special left ideals of the ring of LPDO: 

• chains ([7]) will be finite and different maximal chains for a given L 
they have the same length: if \L) C I\ C Ii C . . . C Ik C |1), \L) C 
J\ C J2 C . . . C J m C |1), then k = m and one can prove a natural 
lattice-theoretic "similarity" of "factors" in both chain. 

• Irreducible LODO will be still irreducible as LPDO. 

• For dim = 2, ord = 2 (that is for operators with two independent 
variables of order two) a LODO is factorizable in this generalized sense 
(i.e. having a nontrivial chain (J7])) iff it is integrable with the Laplace 
cascade method. We describe this cascade method below in subsec- 
tion 13.31 

• Algebraically, the problem is reduced from the ring Q(x, y) [D x , D y ] to 
factorization in rings of formal LODO with noncommutative coefficients 
Q(x, y, D x ) [D y ] and/or Q(x, y, D y ) [D x ] (Ore quotients); in these rings 
all left and right ideals are again principal ideals. 

The details, rather involved, may be found in [51] . This approach neverthe- 
less suffers from the following problems: 

• The definition of divisor ideals given in [51 J is very technical, not intu- 
itive. 

• No algorithms for such generalized factorization is known. 
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A generalization of this ring-theoretic approach to systems of LPDE was pro- 
posed recently by M. Singer. Another ring-theoretic approach was considered 
in [25]. 

In the next subsection we propose a different, much more intuitive defi- 
nition of generalized factorization. 

3.2 General theory of factorization of arbitrary sys- 
tems of LPDE, approach of abelian categories 

Abelian category SCWE of arbitrary systems of LPDE is defined by its 
objects which are simply systems 



S 



L n ui + . . . + L ls U s — 0, a n(<r t Un n 1 

■ u ij c W • • • i ^n) \_ ly x\ t • • • t -^Xn] i 

Lpiiii + • • • + L ps u s = 0, 



Morphism P : S — > Q of two systems is defined as a matrix of differential 
operators 

{Ul = Pu«i + • • • + Pi s « s , 
(9) 
Vrn Pml^l "I - • • • Pms^si 

Pij £ <5(xi, . . . , x n )[D xl , . . . , P Xn ] with the condition that any solution set 
. . . , u s } of the source system (|SJ) is mapped into a subspace of the solution 
space {f i, . . . , u m } of the target system 



Q: 



MuUl + ' ' • + MlmVm = °' My e Q( Xl , . . . , x B ) . . . , D Xn ] 

M q iVi + . . . + M qm V m = 0, 



(10) 

The standard differential Groebner technique (originally developed in the 
beginning of the XX century as the so called Janet-Riquier theory [261 EHl 
HH1I12]) makes this definition constructive: (jHJ) is a morphism mapping (jSJ) to 
ffTUj) iff for any i the equation k M^P^u^ = is reducible to zero modulo 
the equations of the system (jSJ). 

Again, it is easy to see that this category is abelian: for this it is enough to 
check that SCPT>£ is embeddable into the category of (infinite-dimensional) 
vector spaces and linear morphisms and this embedding preserves direct 
sums, products etc. 

It seems natural to refer to Theorem H] to transfer the many properties 
of factorization proved for the category COVO in Section [5] to the case 
of the category SCWE. Unfortunately this is not so simple: the ascend- 
ing chains of monomorphisms are infinite in general: the same example 
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I Ac) C \D 9 ,D™) C ID^D™- 1 ) C ...\D x ,D y ) C |1) makes this obvious. 
The solution to this problem is given by the standard construction of a Serre- 
Grothendieck factorcategory. We refer to [El HU [19] and especially to [T7j 
for a detailed explanation of this important and general construction. One of 
the important steps of this construction is the construction of inverses of mor- 
phisms with "relatively small" kernels; the objects are not formally changed 
in contrast to the ring-theoretic construction of factorrings and factormod- 
ules. In our case we proceed as follows: for a given (say, determined) system 
of LPDE of the form (jSJ) (with s — p) and take the subcategory S n -2 of 
(overdetermined) systems with solution space parameterized by functions of 
at most n — 2 variables. Then the Serre-Grothendieck factorcategory <S/<S n _2 
has finite ascending chains. Another remarkable feature of this factorcate- 
gory is, as we mentioned above, the possibility to consider morphism which 
had kernels defined by systems from S n -2 as invertible morphisms. This may 
lead to a more general theory of Backlund-type transformations (at least for 
the case of linear systems), for example of transformations of Moutard type 

(«). 

Now we can transfer all theoretical results proved in Section [5] to the case 
of the factorcategory S/S n -2 and provide a theoretical foundation for the 
factorization theory of arbitrary linear systems of LPDE. 

The obvious drawback still lies in the absence of algorithms for such a 
generalized factorization. We give an overview of currently known numerous 
partially algorithmic results in the next Sections I3.3H3.51 

3.3 dim = 2, ord = 2: Laplace transformations and Dar- 
boux integrability of nonlinear PDEs 

Here we expose the basics of the classical theory [9j [151, EQ], which is ap- 
plicable to hyperbolic linear partial differential equations of order two with 
two independent variables. For simplicity only the case of an equation with 
straight characteristics will be discussed here: 

Lu = u xy + a(x, y)u x + b(x, y)u y + c(x, y)u = 0. (11) 

The more general case can be found in [2D] El E2]- If one of the Laplace 
invariants of ( iTTjl h = a x + ab — c, k = b y + ab — c vanishes, one can "naively" 
factorize the operator in the l.h.s. of ( |TT1) : k = L = (D y + a) (D x + b); 
h = 0^ L = {D x + b){D y + a). 

If h 7^ 0, k ^ then fill I) is not factorizable in the "naive" sense. In this 
case one can perform one of the two Laplace transformations (not to be mixed 
with Laplace transforms!) which are invertible differential substitutions 
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(isomorphisms in the category of SCWS): 

u = ^(D x + b) tt(i) 

or 

u = - (Dj, + a) w ( _i). 

In fact each of the above substitutions is the inverse of the other up to a 
functional factor. Each of these substitutions produces a new operator of the 
same form (TTTTl but with different coefficients and Laplace invariants. The 
idea of the Laplace cascade method consists in application of these substitu- 
tions a few times, obtaining the (infinite in general) chain 

...<— £(_2) <— £(-i) *— L - *■ £(i) ~^ ^(2) —>•... (12) 

In some cases (namely these cases are considered as integrable in this ap- 
proach) this gives us on some step an operator Lu\ with vanishing hu\ or 
kn\. Then this chain can not be continued further in the respective direction 
and one can find an explicit formula for the complete solution of the trans- 
formed equation; performing the inverse differential substitutions we obtain 
the complete solution of the original equation (with quadratures). 
One of the main results of [5T] are the following Theorems: 

Theorem 5 1 = D x ■ D y — a(x, y)D x — b(x, y)D y — c(x, y) has a nontrivial 
generalized right divisor ideal (so is factorizable in the sense described in 
Section [3|) iff the chain (TTSj) of Laplace transformations is finite at least in 
one direction. 

Theorem 6 1 = D x ■ D y — a(x, y)D x — b(x, y)D y — c{x, y) is a ILCM of two 
generalized right divisor ideals iff the chain fTB) of Laplace transformations 
is finite in both directions. 

This shows the meaning of the generalized definition of [51 J and provides a 
partial algorithm for generalized factorization for equations of the form fill I) . 

Although practically efficient for simple cases, this method has the obvi- 
ous decidability problem: given an operator L, how many steps in the chain 
(112j) should be tried? Currently no stopping criterion is known. As the ex- 
ample (J3J) shows, the number of steps in the chain (equal to n for (J3j) in the 
integrable case c = n(n + 1)) depends on some subtle arithmetic properties 
of the coefficients. 

There exists a remarkable link of the theory of Laplace transformations 
to the theory of integrable nonlinear partial differential equations. This 
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topic was very popular in the XIX century and led to the development of 
integration methods of Lagrange, Monge, Boole and Ampere. G. Darboux 
[TO] generalized the method of Monge (known as the method of intermediate 
integrals) to obtain the most powerful method for exact integration of partial 
differential equations known in the last century. 

Recently in a series of papers [21 HEJ EZ] the Darboux method was cast 
into a more precise and efficient (although not completely algorithmic) form. 
For the case of a single second-order nonlinear PDE of the form 

u xy = F(x,y,u,u x ,u y ) (13) 

the idea consists in linearization: using the substitution u(x,y) — > u(x,y) + 
ev(x, y) and cancelling terms with e n , n > 1, we obtain a LPDE 

v xy = Av x + Bv y + Cv (14) 

with coefficients depending on x, y, u, u x , u y . Equations of the type ( 1141) 
are in fact feasible to the Laplace cascade method, certainly one needs to 
take into consideration the original equation ffTB"]) while performing all the 
computations of the Laplace invariants and Laplace transformations: ffTB"]) 
allows us to express all the mixed derivatives of u via x, y, u and the non- 
mixed u x ... x , u y ...y). The following statement can be found in [20j, recently it 
was rediscovered in [2j HE] : 

Theorem 7 A second order, scalar, hyperbolic partial differential equation 
n~3\) is Darboux integrable if and only if the Laplace sequence IpJ^) for [TJ\ ) 
is finite in both directions. 

In [21 HE] this method was also generalized for the case of a general second- 
order nonlinear PDE 

I/, U, U X , Uy, U XX , U X y, Uyy) 



3.4 dim = 2, ord > 3: Generalized Laplace transforma- 
tions 

In [52] we have proposed a generalization of the Laplace cascade method for 
arbitrary strictly hyperbolic equations with two independent variables of the 
form 

Lu= PiA^y)D x D y u = 0, (15) 

i+j<n 

12 



as well as for n x n first-order linear systems 

n n 

(vi) x = y^Oifc(ar, y)(vk) y + y~]b ik (x, y)v k (16) 
k=i k=i 

with strictly hyperbolic matrix (a^). 

Here we demonstrate this new method on an example of the constant- 
coefficient system 

D x u\ = ui + 2m 2 + m 3 , 

D y u 2 = —6ui + u 2 + 2m 3 , (17) 
(D x + D y )u 3 = 12«i + Qu 2 + m 3 . 

It has the following complete explicit solution: 

ui = 2e^G(x) + e x (3F(y) + F\y)) + exp ^H(x - y), 
u 2 = e y G'(x) + 2e x F'(y) - 2u h 
us = D x u x + 3«i - 2(e y G"(x) + 2e x F / (?/)), 

where F(y), G(x) and if (a; — y) are three arbitrary functions of one variable 
each. 

The solution technology (cf. [52]) for the details) is again a differential 
substitution; in the case of the system (|17|) the transformation is given by: 

Ui = Mi, 

u 2 = u 2 + 2u l , (18) 
us = {{D x + D y )ui - Ui- 2u 2 - 4ui). 

The transformed system has a triangular matrix and is easily integrable: 

D x u 3 = us, 
D y u 2 = 2u 3 + u~2, 
(D x + D v )u\ = us + 2u 2 + Mi. 

Again no stopping criterion for the sequences of generalized Laplace trans- 
formations is known in the general case. For constant coefficient systems an 
alternative technology was proposed by F.Schwarz (private communication, 
2005): transform the system (I17p into a Janet (Grobner) normal form with 
term order: LEX, u 3 > u 2 > Ui,x > y: 

U\ yXX y U\^ XX + U\^ X yy 'iUl^y + 2U\^ X U\^yy + 2U\ ^y U\ = 0, 

u 2 , y + 3u 2 - 2ui jX + 8ui = 0, 

u 2 ,x ~u 2 - \u hxx - \u 1)XV + 3u 1)X + \u Xj y - |«i = 0, 
m 3 + 2m 2 — Ui )X + Ml = 0. 

13 



The first equation factors: 

D 2 x D y -D 2 X + D x D 2 y - 3D x D y + 2D X - D 2 y + 2D V - 1 

= (D x + D y -l)(D y -l)(D x -l). 

So one can find U\ easily and then the other two functions u 2 and are ob- 
tained from the remaining equations of the Janet base producing essentially 
the same solution ( 1T81) . 

Conjecture: For constant-coefficient systems this Grobner basis tech- 
nology is equivalent to the generalized Laplace technology. 

3.5 dim > 3, ord = 2: Dini transformations 

In [TT] another simple generalization of Laplace transformations formally 
applicable to some second-order operators in the space of arbitrary dimension 
was proposed. Namely, suppose that an operator L has its principal symbol 

Sym= a hi2(x)D Xn D Xi2 

h+i2=2 

which factors (as a formal polynomial in formal commutative variables D x .) 
into product of two first-order factors: Sym = XiX 2 (Xj = £\ bij(x)D x . are 
first-order operators) and moreover the complete operator L may be written 
at least in one of the characteristic forms: 

L = {XtX 2 + aiXx + a 2 X 2 + a 3 ) 
L = (X 2 X 1 + aiX 1 + a 2 X 2 + a 3 ), 

where «j = ai(x,y). Since the operators Xj do not necessarily commute we 
have to take into consideration in ( JT9l) and everywhere below the commuta- 
tion law 

[X U X 2 ] = X,X 2 - X 2 X X = P(x, y)X x + Q(x, y)X 2 . (20) 

. This is very restrictive since the two tangent vectors corresponding to 
the first-order operators Xi no longer span the complete tangent space at a 
generic point (xq). ( 1201) is also possible only in the case when these two vec- 
tors give an integrable two-dimensional distribution of the tangent subplanes 
in the sense of Frobenius, i.e. when one can make a change of the independent 
variables (x) such that Xj become parallel to the coordinate plane (xi,x 2 ); 
thus in fact we have an operator L with only D Xl , D X2 in it and we have got 
no really significant generalization of the Laplace method. If one has only 



14 



(TH?]) but f[2"Ul) does not hold one can not perform more that one step in the 
Laplace chain (fT2l and there is no possibility to get an operator with a zero 
Laplace invariant (so naively factorizable and solvable). 

Below we demonstrate on an example, following an approach proposed 
by U. Dini in another paper [12], that one can find a better analogue of 
Laplace transformations for the case when the dimension of the underlying 
space of independent variables is greater than two. Another particular special 
transformation was also proposed in [3], [56]; it is applicable to systems whose 
order coincides with the number of independent variables. The results of [3] , 
[56] lie beyond the scope of this paper. 

Let us take the following equation: 



It has three independent derivatives D X} D y , D z , so the Laplace method is 
not applicable. On the other hand its principal symbol splits into product of 
two first-order factors: £i£ 2 + #£i£3 = 6i (62 + ^63)- This is no longer a typical 
case for hyperbolic operators in dimension 3; we will use this special feature 
introducing two characteristic operators X\ = D x , X 2 = D y + xD z . We have 
again a nontrivial commutator [Ai, A2] = D z = A3. The three operators Aj 
span the complete tangent space in every point (x, y, z). Using them one can 
represent the original second-order operator in one of two partially factorized 
forms: 



Let us use the first one and transform the equation into a system of two 
first-order equations: 



Cross-differentiating the left hand sides of ( l22i) and using the obvious identity 
[X U X 3 ] = [D X ,D Z ] = we get X x X 2 v = D x (D y + xD z )v = X 3 v = D z v or 
= D x (D y + xD z )v-D z v = (D x D y + xD x D z )v = (D y + xD z )D x v = X 2 X x v. 

This is precisely the procedure proposed by Dini in [T2]. Since it re- 
sults now in another second-order equation which is "naively" factorizable 
we easily find its complete solution: 



where <p and ip are two arbitrary functions of two variables each; they give 
the general solutions of the equations X 2 (p = 0, X\ip = 0. 





L — X 2 X\ — X% — X\X 2 — 2 A3. 




(22) 



v = 
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Now we can find u: 



u 



J (vdx+(D y + xD z )vdz) + 0(y), 



where an extra free function 9 of one variable appears as a result of integration 
in (22). 

So we have seen that such Dini transformations (1221) in some cases may 
produce a complete solution in explicit form for a non-trivial three-dimensional 
equation f[2"TT) . This explicit solution can be used to solve initial value prob- 
lems for fl2T|) . 

Dini did not give any general statement on the range of applicability of 
his trick. In [53] we have proved the following 

Theorem 8 Let L = ^2 i+ j +k<2 a^ix, y, z)D x D^D k z have factorizable prin- 
cipal symbol: J2i+j+k=2 a ijk(x, y, z )D x D^ y D k z = S1&2 (mod lower-order terms) 
with generic (non-commuting) first-order LPDO Si, S\. Then there exist 
two Dini transformations Lm, -k(-l) of L. 

Proof. One can represent L in two possible ways: 

L = SiS 2 + f + a(x, y, z) = S 2 Si + U + a(x, y, z) (23) 

with some first-order operators T, U. We will consider the first one obtaining 
a transformation of L into an operator Lrn of similar form. 

In the generic case the operators Si, S2, T span the complete 3-dimensional 
tangent space in a generic point (x,y,z). Precisely this requirement will be 
assumed to hold hereafter; operators L with this property will be called 
generic. 

Let us fix the coefficients in the expansions of the following commutators: 
[S 2 , f] = K(x, y, z)4 + M(x, y, z)S 2 + N(x, y, z)f. (24) 

[Si, S 2 ] = P{x, y, z)Si + Q(x, y, z)S 2 + R{x, y, z)f. (25) 

First we try to represent the operator in a partially factorized form: L = 
(Si + a)(S2 + (3) + V + b(x,y,z) with some indefinite a = a(x,y,z), (3 = 
(3(x, y, z) and V = T — j3S\ — aS 2 , b = a — a[3 — Si(fl). 

Then introducing v = (S 2 + fi)u we get the corresponding first-order 
system: 

|^ (V + b)u = -(Si + a)v. 
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Next we try to eliminate u by cross-differentiating the left hand sides, which 
gives 

[(V + b),(S 2 + (3)}u = (S 2 + /3)(S 1 + a)v + (V + b)v. (27) 

If one wants u to disappear from this new equation one should find out when 
[{V + b), (S 2 + (3)}u can be transformed into an expression involving only v, 
i.e. when this commutator is a linear combination of just two expressions 
(S 2 + (3) and (V + b): 

[(V + b), (S 2 + (3)] = ti(x, y, z){S 2 + (3) + v(x, y, z)(V + b). (28) 

This is possible to achieve choosing the free functions a(x,y,z), (3{x,y,z) 
appropriately. In fact, expanding the left and right hand sides in (f28l) in 
the local basis of the initial fixed operators Si, S 2 , T and the zeroth-order 
operator 1 and collecting the coefficients of this expansion, one gets the 
following system for the unknown functions a, (3, /1, v: 

K + (3P-S 2 {(3) = up, 
M — S 2 (a) + j3Q = va — \i, 
N + pR= -v, 

(313^(3) - f{(3) + S 2 {a) - (3S 2 (a) - S 2 (S x (p)) = -v{a - a/3 - - ///3. 

After elimination of v from its first and third equations we get a first-order 
non-linear partial differential equation for f3: 

S 2 {(3) = (3 2 R+{N + P)(3 + K. (29) 

This Riccati-like equation may be transformed into a second-order linear 
PDE via the standard substitution (3 = S 2 {'i)/^(- Taking any non-zero solu- 
tion j3 of this equation and substituting fi = ua + S 2 (a) —(3Q — M (taken from 
the second equation of the system) into the fourth equation of the system 
we obtain a first-order linear partial differential equation for a with the first- 
order term (3 S 2 (a). Any solution of this equation will give the necessary value 
of a. Now we can substitute [(V + b), (S 2 + (3)\u = n(S 2 + (3)u + v{V + b)u = 
fiv — v{S\ + a)v into the left hand side of ( 1271) obtaining the transformed 
equation L^v = 0. 

If we would start the same procedure using the second partial factorization 
in (T2U]) we would find the other transformed equation L(_!)U> = 0. □ 

4 Other results and conjectures 

The theory of integration of linear and nonlinear partial differential equations 
was among the most popular topics in the XIX century. Enormous amount 



17 



of papers were devoted for example to transformations of equations to an 
integrable form. In particular the papers EH EE] were devoted to a more 
general Laplace type transformations. Some of these results were obtained in 
the framework of the classical differential geometry; cf. [TH [27] for a modern 
exposition of those results. 

In addition to the problems studied above one should mention a class 
of overdetermined systems of linear partial differential equations with finite- 
dimensional solution space studied in [58l [32] ■ There an algorithm for fac- 
torization of such systems was proposed. 

Another popular in the past decade topic was the theory of "naive" fac- 
torization, i.e. representation of a given LPDO as a product of lower-order 
LPDO: in [21] an algorithm for such factorization was proposed for the case 
of operators with symbol representable as a product of two coprime polyno- 
mials. This result was developed further in [15] . 

From the theory of Laplace and Dini transformations the following con- 
jectures seem to be natural: 

• If a LPDO is factorizable in the generalized sense, then its principal 
symbol is factorizable as a multivariate commutative polynomial. 

• If a LPDO of order n has a complete solution in a quadrature-free form 
(TJP then its symbol splits into n linear factors. 
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